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Resonance in an open quantum dot system with a Coulomb interaction: 

a Bethe-ansatz approach 
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An open quantum system consisting of a quantum dot with a Coulomb interaction and 
two leads without interactions is studied. The many-body scattering states are constructed 
with the Bethe-ansatz approach. The expectation value of the electric current is exactly 
calculated for the scattering states to observe resonance peaks due to many-body scattering. 
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The purpose of this letter is to observe resonance in an open quantum system with a 
Coulomb interaction. The system that we study is the two-lead interacting resonant-level 
model (IRLM) , which consists of two leads of non- interacting electrons that interact with an 
electron on a quantum dot in between the two leads. We obtain iV-electron scattering states for 
arbitrary N, generalizing the Bethe-ansatz approach to open systems. By using the scattering 
states, we exactly calculate the quantum-mechanical expectation value of the electric current 
through the quantum dot, thereby observing resonance peaks. Some of the resonance peaks 
appear only when the interaction exists; they reflect the effect of many-body scattering. 

The resonance of many-body scattering that we observe in the quantum-mechanical ex- 
pectation value has not been found in previous works with the Bethe ansatz. The Bethe-ansatz 
approach has provided a nonperturbative method of studying equilibrium states of interact- 
ing quantum systems including the Kondo problem. 1-5 The approach is now used to discuss 
transport properties of mesoscopic systems. Konik et al. 6,7 studied transport properties of the 
Anderson model in the thermodynamic limit of a closed system with periodic boundary con- 
ditions. Our scattering states, in contrast, appear only in open systems; they are constructed 
without imposing periodic boundary conditions. By extending the Bethe-ansatz approach, 
Mehta and Andrei 8 studied the two-lead IRLM as an open system to obtain iV-electron scat- 
tering states giving nonequilibrium steady states in the limit TV — > oo. In their study, however, 
the quantum-mechanical expectation value of the current does not depend on the interaction; 
the effect of the interaction appears only in the statistical-mechanical expectation value as 
modification of the Fermi distribution in the leads. Thus our results are different from the 
previous ones. 
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There has recently been a great deal of interest in mesoscopic systems with interacting 
electrons. Experiments suggest that interactions are essential in understanding their transport 
properties. 9-13 The perturbation theory tells us that the effect of interactions is observed as 
resonance peaks of the electrical conductance. 14 ' 15 For non-interacting open quantum systems, 
the relation between quantum mechanical scattering states and nonequilibrium steady states 
is well investigated. 16 ' 17 However, the relation in interacting open quantum systems has not 
been clarified, excepting Schiller and Hershfield's result 18 at a special point of the interaction 
parameter where an interacting system is mapped to a non-interacting one. The present 
study gives a steady step toward an exact analysis of interacting open quantum systems out 
of equilibrium. 

The Hamiltonian of the two-lead IRLM is given by 

H = E ( pM^ z) \^ z) + ^(4(0)^+^^(0))) 

+ e d dU + U 4(0)ci(0)dU, (1) 

i=l,2 

where i(> 0) is the transfer integral between each lead and the quantum dot, is the gate 
energy of the dot and U(> 0) expresses the Coulomb repulsion. The dispersion relation in 
the leads is linearized in the vicinity of the Fermi energy to be E = k, under the assumption 
that t, Ed and U are small compared with the Fermi energy. 1 ' 3 ' 5 We stress that we treat the 
system as an open system in the limit L — > oo. The one-lead IRLM with periodic boundary 
conditions was studied with the Bethe ansatz. 3 Our purpose is to investigate, for scattering 
states, the electric current through the quantum dot, 

I:=N 2 -N 1 = ± ^(-)*(4(0)d-^(0)). (2) 



i=l,2 



We derive the Schrodinger equations for the system. After the transformation c 1 / 2 (-z) = 
(c e (z) ± c (z)) /\/2, the Hamiltonian (1) is decomposed into the even and odd parts. 8 Due to 
the relations [H,N e + Nd] = [H,N a ] = for the number operators N e / Q = f dzc . (z)c e / (z) 
and Nd = d^d, the set (N e + Nd, N a ) gives a good quantum number. The iV-electron state in 
the subspace with N a = n, (0 < n < N) is generally expressed in the form 

l^> = (J dz g {n) \z)c\{z{) ■ ■ ■ cl(z N -n)cl(z N -n+l) ■ ■ ■ c\{z N ) 

+ fdze^{z)cl( Zl ) ■ ■ ■ cUzN-n^dtcUzN-n) • • • c\{z N ^)) |0), (3) 

where g {n) (z) = ( Zl , . . . , z N ) for < n ^ N and e {n \z) = e^ n \ Zl , . . . , z N -i) for < n < 
N — 1 are functions to be determined. We also set (z) = for convenience. The eigenvalue 
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problem H\ip) = E\ip) is cast into a set of the Schrddinger equations 

N „ N—n 

(N - „)(]T - — -£)<,<»>(,) + ( Y, (-) N -'-S(z,)e"»(. . . , z,. u z, +u . ..)=0, 



1=1 

JV-1 



( E (TdI + U ^ Zi) ) +e d -EY n \z)+t{N-n)g { - n \. . . , z N - n - U 0, z N - n , . . .) = 0, (4) 
«=i 

where O^n^iV— 1. In what follows, we use the variables and j/j to express the 
coordinates of the leads e and o, respectively: g^ n \z) = g( n \x\, . . . , xjv-n, 2/i> • • • , Vn) and 
e ( n )( z ) = e( n )(xi, . . . , xjv-n-ii Z/i) ■ ■ ■ j Vn)- The set of eigenfunctions in the one-electron sector 
with E = k is given by 

5 (0) (xi) = 0fc(si) := -^^(fl(-xi) + e^(*i)), 



e (0) = e fe : = 



k-e d ' 



5 (1) (yi) = ^(yi):= " 



1 + e i5fe 

with the phase shift := —2 arctan (t 2 /2(k — e^)) of one-body scattering at xi = in the lead 
e and the step function 6{z). Note that the eigenfunction g k {x\) is discontinuous at x\ = 0. 

We construct an iV-electron eigenstate with the Bethe ansatz. It is different from the one 
obtained by Mehta and Andrei. 8 To demonstrate the difference, we first consider the case 
N = 2. The set of two-electron eigenfunctions with the energy eigenvalue E = k\ + k 2 is 
assumed to be 

2 4ifc 2 ( x i' x 2) = Z fclfc2 (xi - x 2 )g kl {xi)g k2 {x 2 ) - Z klk2 (x 2 - x 1 )g k , 2 (x 1 )g kl (x 2 ), 
4il 2 ( Xl ) = Z kik 2 ( x i)9k 1 (x 1 )e k2 -Z klk2 (-x 1 )g k2 (x 1 )e kl , 
flSafal'fl) = X ( x l -yi)5fci(«l)^fca(j/l)> 

4i 2 (yi) = ^(-yi) e fci /i fc 2 (yi)> 

2 9 k % 2 (Vi,y2) = hk 1 (yi)h k2 (y 2 )-h k2 {y 1 )h kl (y 2 ), (5) 
where the amplitudes Z klk2 (z) and are defined by 

£fcife,(*) := e-^ fc i fc 20(-z) +e^ k ^e(z), 
X(z) := e^0(-2) +e~^0( z ) 



with the phase shifts p klk2 := 2 arctan ( — an d ?? := 2 arctan(— J7/2) of two-body 

scattering. The eigenfunctions 5^ (zi, £2) and e^ 2 (^i) are the same Bethe eigenfunctions 
as those assumed in the one-lead IRLM, 3 although we do not impose periodic boundary 
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conditions. The eigenfunctions g^ k2 {x\,yi) and (yi) are obtained with separation of 
variables; if we set yj^On, yi) = X(xi - yi)g { ^ k2 (xi,yi) and ejJJ^yi) = X^y^e^^i), 
Eqs. (4) are decoupled into even and odd parts, and the eigenfunctions 9^k 2 {xi,yi) and 
(j/i) are given by the product of eigenfunctions of the even and the odd parts. The 

(2) 

eigenfunctions g klk2 (yi, 2/2) should be a free fermion eigenfunction because of Eq. (4) for 
g( N \z). The phase shifts in our solution are different for each two-body scattering in the lead 
e, in the lead o and between the leads; this gives the resonance of many-body scattering in 
the current expectation value, as we shall see below. In Mehta and Andrei's solution, 8 on 
the other hand, the same phase shift <p klk2 °f two-body scattering was adopted for all the 
two-electron eigenfunctions and e^ n \z). 

By exchanging k\ and k% in both eigenfunctions g k ^ k {x±,yi) and e^ k {vi ) > we have another 
set of eigenfunctions g^ kl (xi,yi) and ejj.*^. (yi) with the same eigenvalue E = k± + &2. In the 
limit t,U -> 0, the set {yj^ (xi, x 2 ), yj^fci, yi), yj^fci (^1, yi), yj£i 2 (yi, y2)} reproduces a 
complete orthogonal system of two free fermions in the two leads, while Mehta and Andrei's 
solution 8 does not. In this sense, our solution (5) is more plausible than theirs. 

In a way similar to the case N = 2, we obtain a set of TV-electron eigenfunctions with the 
energy eigenvalue E = YliLi i n the form 

Sfc^fof) = fjV-nVn' 2 sgn(PQ)4" ) (^;y) 5fc Pi (a*) JJ (%■), 

Qes„ 

j N—n—1 n 

e { k \x;y) = (N _ n _ 1)w ^2 s MPQ)z ( k n J(x;y) JJ gk Pi (xi)e kPN _ n Y[h kN _ n+Qj ( yj ), (6) 

1 '' '-P6SjV-n 1=1 j = l 

Qes„ 

where & m is the symmetric group acting on the set {1,2,..., m} and 

Z k n J( x ;y)= II Z k Pi k P .(xi -Xj) YI X(xi-yj), 

1 ^ j ^ n 

i£ ) (x;y) = z£)(x;y)| x ^ Ml . 

Note that the amplitude Z k n J(x; y) and zj^(x; y) are given by the product of the amplitudes 
of the two-electron eigenfunctions. The iV-electron eigenfunctions (6) are indexed by a set 
of momenta k = {fei, . . . , fejy}. We denote by \k;n) the eigenstate obtained by putting the 
eigenfunctions g k n \x;y) and e k n \x;y) into (3), and call it a Bethe eigenstate. 

We show that, for a fixed set of N momenta k = {k±, . . . , kw}, there exist 2 N degenerate 
Bethe eigenstates with the energy eigenvalue E = YliLi For a fixed n, we consider jvC„ 
ways of dividing the set k into two subsets wherein the first subset contains N — n elements 
and the second subset contains n elements. It is convenient to index each way of dividing 
by kn = {kn 1} . . . , kn N } with an element R of the symmetric group &n satisfying R\ < 
R2 < ••• < Rn-u and Rn-h+i < Rn-u+2 < ••• < Rn- The element R is an element of 
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&N /(&N-n x ®n)) where &N-n is the symmetric group acting on {1,2,- •• , N — n} and & n 
that acting on {N - n + 1, N - n + 2, . . . , N}. For ^ n ^ N and i? G 6Ar/(6Ar„„ x 6„), 
all the Bethe eigenstates |&R;n) with a set of momenta have the same energy eigenvalue 
E = YliLi ki- In the limit L — > oo, the Bethe eigenstates |/c.R;n) satisfy the relation 

N 2 

(k R ; n\k s ;m) = 5 nm 5 RS L N ]J — — + O^' 1 ) 

■ i J- ~r cos Ok 

for generic values of {ki}. Hence the normalized Bethe eigenstates are orthogonal in the limit 
L — > oo. As a result, the total degree of degeneracy of the energy eigenvalue E = YliLi ki is 

We obtain a general iV-electron eigenstate by taking a linear combination of the 2^ de- 
generate Bethe eigenstates \kR-,n) in the form 

N 

llcl^^sgn^l^;,}, (7) 

n=0 R 

where the sum on R runs over elements in ©at/ (&n-u x6 n ). The square norm of the eigenstate 
\k) is readily calculated from (k\k) = Y^h=()Y^r\^r\ 2 {kR;n\kR]n) . 

The expectation value (I) = (k\I\k) / (k\k) of the current operator / in (2) for each eigen- 
state \k) in (7) is exactly given by 

N 

V> 2 N 'iL ^ (N-n)\n\ Z^ in H A P A P e > 

n=i v 7 Pe6iv 

x (ll cos 2 J v COS 2J C ° S 2 efc ^-„ + i+°( L )• ( 8 ) 

1=1 

Here, by using the fact that any element P G &n is decomposed as P = RQ with a unique 
element R G ©^/(©jv-n x ©n) an d a unique element Q G ©v-n x ©n 5 we set ylp^ = yl^ 
for every P with the same R. 

The result in (8) shows that we need to consider the linear combination (7) of the degen- 
erate Bethe states in order to obtain a non-zero current expectation value. Indeed, a specific 
Bethe eigenstate \kn;n) gives the expectation value [kR] n\I\ku; n) = 0. We stress that we do 
not impose periodic boundary conditions to the eigenfunctions g^\x;y) and e^ l \x;y) in (6); 
with the periodic boundary conditions, 6 ' 7 the eigenstates \kn;n) with different n or different 
R would not be degenerate, and hence the current expectation value would be zero. 

By expressing the eigenstate \k) in terms of the leads 1 and 2, the eigenfunction describing 
N — n electrons in the lead 1 and n electrons in the lead 2 is given by 

(iV-n)!n! F^ n \z) := (c 2 (z N ) • • • C2(z N - n +i)ci(z N - n ) • • • ci(^i)|fe) 

1 N 

= _L_ ^ y^^-)i{ R i\ N - n < R N-^+l,-,RN} sgn ^ 

2 2 m= o R 
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x {c (zr n ) ■ ■ ■ c (zR N _ m+1 )c e (z RN _J ■ ■ ■ c e (z Rl )\k), 

where §A stands for the number of elements in the set A. We consider the behavior of (z) 
in the region z\ < Z2 < • • • < zn < 0. The eigenfunction F^ a \z) is a complicated linear 
combination of plane waves e 1 ^ iA '' p » Zi for P £ Among them, we call the plain wave 
e iJ2i k i z i an "incoming wave". The terms with the incoming wave are summarized as 

e i Ei k i z i y^^_y{R i \N-n<R N - m+1 ,...,R N }^(m) 
m,R 

where 

~Am) Jm) -U^i<3iN^<Pk R k R +ri'EisiN^ m sgn(R i -R j )) 

A R ' = A R 'e 3 N~m<j 

We define the scattering states \k)^\ (£ = 0, 1, . . . , N) by taking the coefficients {^4^} of 
the eigenstate \k) as 

J[M = (_y{Ri\ N -e<RN-n+i,-,RN} ; (g) 

for ki > k,2 > • • • > k]y- In the scattering state \k)( e \ the incoming wave e 1 ^*^!^ exists only 
in the eigenfunction Fjf\z), which describes N — £ electrons in the lead 1 and £ electrons in 
the lead 2. 

The scattering states |fe)^ are different from those of the standard one-body scattering 
theory in quantum mechanics. If we were solving the one-body scattering problem, the scat- 
tering state would be obtained from the condition that an electron comes only from the lead 
1 or the lead 2. However, the eigenstate \k) in (7) does not give such scattering state. In fact, 
the scattering state |fe)^ extends to all parts of the two leads for U > 0. In other words, it is 
impossible to judge whether each electron comes from the lead 1 or the lead 2 for U > 0, which 
is not strange since we assume the same Fermi energy for both leads. In the limit U — > 0, our 
scattering state \k)^ is reduced to the standard one-body scattering state. 
By applying (9) to the expectation value (8), we have 

N 

t \ - n 



n=l ( ' PeSjv 

j n-l N 
-{^sgn{Pi-P n )ip kPikpn +^sgn(Pj-P n )?7- S k 
i=i i=i 

/TT ^kp+V ^f 7]\N-n 6 k 

U1 C ° S 2 Jl COS 2/ cos l" e ^ + °^ )■ 



x 

i=i 



Short calculations reveal that every term in (1)^ contains the product of the factors sin 5 k , 1 + 
cos S k , sm((p kikj + rf) and 1 + cos(<p kikj + 77), which are rational functions of ki, t, td and U. The 
factors have poles at = k{ ± it 2 /2, (ki + kj)/2 =F i{ki — kj)U/4 in the complex plane of en- 
Figure 1 shows the current expectation value (1)^ as a function of the gate energy for the 
scattering states \k) (t > indexed by (N,£) = (2,0) and (3, 1). 
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Fig. 1. The current expectation value (J)W for the scattering states. We fixed t = 1. 



We find resonance peaks in the vicinity of = (ki + kj)/2, which correspond to many- 
body scattering; they appear only for U > 0. As is stressed above, the resonance of many-body 
scattering is originated from the phase shifts which are different for each two-body scattering 
in the lead e, in the lead o and between the two leads. We also find resonance peaks in the 
vicinity of = ki, which correspond to one-body scattering at the quantum dot and are 
reduced to Lorentzian peaks in the limit U — > 0. The resonance peaks in the vicinity of 
e d = (ki + kj)/2 were not present in Mehta and Andrei's result; 8 their results are equal to 
the limit U — > of our result. This is because the interaction effect would be canceled in 
the current expectation value {1} if we adopted the same phase shifts for all the two-body 
scattering in the lead e, in the lead o and between the two leads. Our choice (5) of the phase 
shifts of two-body scattering is more plausible in the context of eigenstates as mentioned 
above. It would be interesting to discuss how the resonance of many-body scattering affects 
the transport properties of the interacting open quantum system out of equilibrium. 
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